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ANALISE MATEMÁTICA IV 

FICHA SUPLEMENTAR 3 

EQUAÇÕES DIFERENCIAIS ESCALARES DE PRIMEIRA ORDEM 

Equações Lineares Homogéneas 




Resolução: Para y ^ 0, 



y e 



«=> f±dy = -J édt + c 

<í=> In \y\ = —ê + c 

<^=> \y(t)\=ke~ et ondek>0 

<í=> y(t) = ke~ e onde k ^ 

Quando y = 0, encontra-se que a função y{t) = , Vi e R, também é solução. 
Solução geral: 

y(t) = ke~ é com k <E R . 

Intervalo de definição: R. 
Verificação: 



logo 



^L = ±(ke- et ) = ki-e^e^ 
dt dt 



ok! 



D 



Comentário: Esta EDO tem uma família infinita de soluções, família essa parametrizada 
por k 6 R, ou seja, para cada k real tem-se uma solução. Cada solução está definida em 
toda a recta real. (} 



Determine a solução do seguinte problema de valor inicial 



(2) 



§ + vIT^y = o 



dt 



2/(0) = \/5 . 



Resolução: 

Resolução da EDO : Para y ^ 0, 



y 



hi\y\ = - f Vi + t 2 dt 
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c— x def 



Cálculo de uma primitiva de Vi + t 2 : Fazendo a substituição t = - — j — = sinhx, para a qual 
= cosh x, obtém-se 



d.r 



J Vi + í 2 dt = J" Vi + sinh a; cosh a; da; 



cosh a; 

,2 



= J cosh xdx a qual se pode primitivar por partes, ficando 
= sinh x cosh x — j sinh x dx 

cosh 2 x— 1 



Das duas últimas linhas obtém-se 



2 / cosh xdx = sinh a; cosh x + x + c 



pelo que 



V 1 + t 2 dt = — (sinh x cosh x + x) + c . 



Fa/ta escrever esta primitiva em termos de t: 

t= eX -f* ^^ (e x ) 2 - 2te x - 1 = 

^^ e* = t + Vl + í 2 



onc/e se aplicou a fórmula resolvente para a equação quadrática (t — Vi + t 2 não pode ser solução 
porque é sempre negativo). Logo, 

x = ln(í + Vl + í 2 ) . 
Por outro lado, 

sinh x = t e cosh x = vl + í 2 . 

Finalmente, a primitiva em questão fica 

jVTTt^dt = i(ívm 2 +in(í+vm 2 )) + c 



= iíVl + í 2 + ln y/t + Vi + í 2 " + c . 
Voltando à EDO, 

]n\y\ = - f \/l + t 2 dt= -\t\/i + í 2 - ln Vi + Vi + í 2 + c 

<^> |y(í) | = fc exp(- \t\/\ + í 2 - ln Vi + Vi + í 2 ) onde k > 

^> y(t) = k exp(- iíVl + í 2 - ln Vi + Vi + í 2 ) onde k ^ Q 

^^ y (t) = k-r-± e -^^ /I+15 " ondek^O 

Quando y = 0, encontra-se que a função y(t) = , Vi e IR, também é solução. 
Solução geral: 

y(t) = k ■ = = e~^^ /I+1T com fc e IR . 

Ví + vTTT 2 

Intervalo de definição: IR. 
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Verificação: 



4í - i [k exp(- ±í VT+t 5 - In Vi + Vi + í 2 ) 
«(*) • i (-htVTT^ - ln Vi + Vi + í 2 ) 



= „(t). -iVTT^-i 



i í 2 



/i+t-= 



2 v / l+t 2 ~ 2(t+VT+F) 

? ,m . r_I,/TT72 _ í 2 (t+VT+t T )+^T+F+A 

/ \ 



V 



i./TTT? (i+t 2 )(t+>/í+P) 
'2 Vi + I 2yi+F(í+yT+^) 



v / I+F / 



= -VI+í 2 ^) 

Condição inicial : V5 = y(0) = &e° 
Solução do problema de valor inicial: 



-ok! 

k = V5. 



y(t) = 



t+ vr+í 2 



e ^ tVI + ¥ onde k e 



D 



Comentário: Csíe problema de valor inicial tem uma única solução definida para todo o 
t real. <> 



Equações Lineares 




Resolução: O factor de integração desta EDO linear é e*. Multiplicando por é tem-se 



+ dy + e 



í (e ' s) = rb 



e integrando entre 2 et, tem-se 



éy{t) - e 2 y{2) = J ^-^d, 



1 + s 2 



y(t) = e-* (te 2 + J - 
y(t) = Se 2 "* + e _i 



t e .s 






:ds 



1 + S 2 



ds 



Solução: 



Intervalo de definição: 



y(í) = 3e 2 - í + e- t (£ T 



t e s 



:ds 
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Verificação: A condição inicial é verificada 



y{2) = 3e° + e°(í j^dsj =3 + = 3 



e a equação também: 



dy_ 

dt 



-3e 2 -* - e"* 



-3e 2 -* - e"* 



t n s 



f , -t d 
ds + e — 



1 + s 2 



t &s 



dt \Jo 1 + s 2 



ds 



= "!/(*) + 



1 + í 2 



^-cís + e 

1 + s 2 



o/c/ 



1 + í 2 



D 



Comentário: Como nem sempre é possível primitivar uma função em termos de funções 
elementares, pode acontecer (como no exercício anterior) que a solução de uma equação 
diferencial só possa ser apresentada em termos de um integral indefinido. <0> 



(4) 



Resolva o seguinte problema de valor inicial 




^L + h(t)y = t , y(-l) = 


= 2 


onde h(t) é a função definida por 




, ,,. f se t < 
h ® = \t seí>0 





Resolução: O factor de integração para esta EDO linear é eJ 'H í ) dí . Uma primitiva de 
h(t) pode obter-se por meio de um integral indefinido, por exemplo: 



H(t) = j h{s)ds = | ? É 2 



se t < 
set > 



Multiplicando pelo factor de integração tem-se 



e m)^ + h{t)e H(t )y = te H(t) 



d 

H 



(e H ^y) = te 



H(t) 
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e integrando entre —1 et, tem-se 

e m y{t) _ e *(-i) y (_i) = /"* se^Wds 

2e -H(t) + e -H(t) ^ se o ds + jt ge l s 2 ds 
2e -H(t) + e -H(t) j-t i se o ds 



!/(*) = 



sei > O 
sei < O 



*(*) = 



*(*) 



y(t) = 



1+2 

26-2* 



e 2' 



-è + 



is 2 

e2 s 



sei > O 



ze + 2 2 



26-2* 



e 2 e2 



ií 2 /Mt 2 



z ^ 2 2 



l„-£i 



sei < O 
se í > O 

sei < O 



e"2 r +1 sei > O 



sei < O 



Solução: 



y(t) = 



\e~^ +1 sei > O 



sei< O 



Intervalo de definição: R 
Verificação: Primeiro nota-se que 

, s 1 3 
^" 1 )=2 + 2= 2 
/4/em d/sso, tem-se, para t ^ O 

-\te~^ = -íy(í) + í = -h(t)y(t) +t se t > O 
í= — Oí/(í) + * = -h(t)y(t) + t 



dy_ 

dt 



set < O 
£m í = O ambas as derivadas laterais de y(t) são O e portanto a equação 



^(0) + h(0)y(0) = 



é satisfeita. 



D 



Comentário: Apesar da expressão da solução do problema de valor inicial ter de ser dada 
por ramos, a função y(t) é de classe C 1 . Isto deve-se ao facto de H(t) ser de classe C 1 
em R (é uma primitiva da função contínua h(i)) e à expressão integral para a solução y(t). 

O 
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(5) 



Um cardume de salmões vive tranquilamente numa zona costeira. A taxa de nata- 
lidade do cardume é de 2 por cento por dia e a taxa de mortalidade de 1 por cento 
por dia. Em t=0 o cardume tem 1000 salmões e nesse instante chega à zona um 
tubarão que se dedica a consumir 15 salmões por dia. Quanto tempo demora o 
tubarão a extinguir o cardume? 



Resolução: Seja y(t) a população de salmões no dia t. A evolução normal da população 
seria aumentar à taxa de 2 — 1 = 1 por cento por dia. Isto é a população deveria satisfazer 
a equação diferencial: 

d -l = 0.01, 
dt 

No entanto, a partir do momento em que o tubarão chega, morrem mais 15 salmões por 
dia e portanto tem-se 

^ = 0.01i/ -15 
dt 

Uma vez que no instante t = há 1000 salmões, para responder à questão do enunciado 
tem de se resolver para t > o problema de valor inicial: 



^ = 0.01y-15 

dt 



?/(0) = 1000 



O factor de integração para esta EDO linear é e 001í . Multiplicando pelo factor de inte- 
gração tem-se 

| (c-°- 01i „(t)) = -15 e -°- 01í 

e integrando entre et tem-se 



-O.Olí 



y(t) - 1000 = 



-15e-°- 01s d S 



y(t) = e°- 01í ( 



1000 



15e-°- 01a d S 



y(t) = lOOOe 



O.Olí 



15e 



o 

O.Olí 



-O.Ols 



0.01 



y(t) = 1000e°- olí - 1500(e a01í - 1) 
y(t) = 1500 - 500e°- olí 



Solução: 



Intervalo de definição: 



y(t) = 1500 - 500e 



O.Olí 



Note-se que apesar de a solução da equação estar definida para íel, ela só tem significado 
físico para t > e y(t) > 0. 
Verificação: Primeiro nota-se que 

y(0) = 1500 - 500 = 1000 



Além disso, tem-se 

= -5e auit = 0.01(2/ - 1500) = 0.01y(í) - 15 



dlJ r- O.Olí 



dt 



ok\ 
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Resposta ao problema: O cardume estará extinto quando y{t) = O 

1500 - 500e a01í = O 
e o.oit = 3 

t = 100 In 3 
í = 109.8 

Conclui-se que o tubarão leva aproximadamente 110 dias a extinguir o cardume de salmões. 

D 



(6) 



Um resíduo industrial é despejado num tanque cheio com 1000 litros de água a uma 
taxa de 1 litro por minuto. A mistura bem homogénea é despejada à mesma taxa. 

(a) Determine a concentração de resíduos no tanque no instante t. 

(b) Quanto tempo leva esta concentração a atingir os 20 por cento? 



Resolução: 

(a) Seja y(t) a quantidade de resíduo no tanque no minuto t em litros. Se o tanque tivesse 
capacidade ilimitada e não houvesse escoamento, a quantidade de resíduo aumentaria 
de acordo com a lei 

dy = 1 

dt 
Como se está a admitir que o resíduo se mistura imediatamente e que o tanque se 
escoa à mesma taxa em que é cheio, a taxa de aumento de resíduo será inferior 
porque uma parte do resíduo se escoará imediatamente. A questão é: que parte? 
Uma vez que a mistura é homogénea, a proporção de resíduo no litro que se escoa em 
cada minuto será idêntica à proporção de resíduo no tanque. Esta proporção é j^ 
portanto conclui-se que a quantidade de resíduo que se escoa por minuto é 0.001y(í) 
litros. Donde y(t) obedece à lei: 

% = 1 - O.OOly 
dt 

No instante t = não há qualquer resíduo no tanque, portanto 

2/(0) = 

Para resolver esta EDO linear multiplica-se pelo factor de integração e 0001í e obtem- 
se: 

j t {e omit y(t)) = e omu 
e integrando entre et tem-se: 

e 0.001t y{t) _ Q= f e 0MU ds 

Jo 
't 

-0.001Í / / „0.001s. 



^ y(t) = e- auuit I / e u - UUi - 5 ds 

^ y{t) = 1000(1 - e -°- 001í ) 

Solução : A concentração c(t) de resíduo é igual à quantidade de resíduo a dividir pela 
capacidade do tanque. Conclui-se que 

c (t) = i _ e - omit s 

Intervalo de definição : R (mas a solução só tem significado para t > 0) 
Verificação: Primeiro nota-se que 

2/(0) = 1000(1 - 1) = 
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Além disso, tem-se 

"" = e -°- 001í = -0.001(y - 1000) = -0.001y(í) + 1 -ok\ 



dt 
(b) A concentração atinge 20 por cento quando a quantidade de resíduo for 200 litros. 



1000(1 - e - u - UUit ) = 200 

e -o.ooit = ± _ 02 

t = -1000 In 0.8 
t = 223.1 

Portanto a concentração leva aproximadamente 3h e 43m a atingir os 20 por cento. 

D 



(7) 



Chama-se 


equação de Bernoulli a uma 


equação 


diferencial da 


fo 


rma 






dy 
dt 


a(t)y = 


b(t)y n 






onde 


n > 


1 é um número natural. 










(a) 


Determine a solução geral desta eqi 

Sugestão: Divida ambos os 


ação. 

lados por 


y n e transform 


5 a 


equação 


(b) 


Ache 


numa equação linear fazendo a 

a solução geral da seguinte eq 


mudança 

uação de 


de variável u = 

Bernoulli 


y L ~ 


-71 






dy . 3 
— + ysmt + y 


sin 2í = 







Equações Separáveis 




Resolução: Para í 7^ 0, 



,lt - 7 — I^ dt = Ij dt + c 



=^> J dy = In \t\ + c 



^^ y(t) = In \t\ + c 

Solução geral: 

y(t) = In |í| + c com c £ IR . 

Intervalo de definição : ] — 00, 0[ otv ]0, +oo[. 

Verificação: 

dy d ,, . . 1 , , 

— = —(In í + c = - para í ^ 0) - o/c/ 

dt dV M ; í VP ^ ; 



D 



Comentário: ,4 equação dada tem duas famílias infinitas de soluções, cada família 
parametrizada por cei. As soluções de uma família estão definidas no semi-eixo aberto 
] — 00, 0[ e as soluções da outra família no semi-eixo aberto ]0, +00 [. <0> 
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Resolução: Factorizando o segundo membro, e dividindo por 1 + y 2 (que nunca é zero), 

d J = l-t + y 2 -ty 2 ^ % 



f = l-t + y 2 -ty 2 ^ £k = (i_ f )(i + 1/ 2) 



y 



= 1 -t 



l+y 2 

Jj^?dy = f(l-t)dt + c 



arctan y = t — \t 2 + c 



para t-h 2 + c^^ + kn,keZ. 



y(t) = tan(í - \t 2 + c) 

±í 2 + c^§ ' 

Solução geral: 

y(t) = tan(í - -í 2 + c) com c G R . 

Intervalo de definição: Os intervalos de definição possíveis são os intervalos máximos con- 
tidos no conjunto 



1 7T 

\ {í G R : t - -t 2 + c + - + kn , k G Z} 



Verificação: 



f = |[tan(í -ií 2 + c)] 

= (l-í)[l + tan 2 (í-ií 2 + c)] 

l-t + y 2 -ty 2 = (l-t)(l + y 2 ) 

= (1 - í) [1 + tan 2 (í - \t 2 + c)] - ok! 



D 



Comentário: Os possíveis intervalos de definição da solução foram apenas apresentados 
implicitamente como intervalos abertos de comprimento máximo contidos no conjunto 

1 7T 

R \ {í G R : t - -t 2 + c + - + kn , k G Z} . 

Por outras palavras, um intervalo de definição concreto terá por extremos duas soluções 
consecutivas da equação quadrática 

In 7T 

-t 2 -t-c+ - + klT = 

com k a variar em Z. Se fosse dada uma condição incial, escolher-se-ia a constante real c 
e determinar-se-ia o intervalo que contivesse o instante inicial. 

Por exemplo, se fosse dada a condição inicial y(0) = 0, a constante c teria que satisfazer 
= y(0) = tanc. Escolhendo c = 0, o intervalo de definição da solução deste problema 
de valor inicial seria o intervalo máximo contido no conjunto 

1 7T 

R \ {í G R : t - -t 2 + - + kn , k G Z} 



10 
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que contém o instante t = 0. Resolvendo a equação quadrática 

t 2 - 2í + 7T + 2kTT = ^> t = 1 ± Vi - 7T - 2/C7T 

("com fc g TL), obtém-se que as soluções mais próximas det = são 

e 



t = 1- VT 



Í = 1 + V1 



7T . 



Sendo estes os dois instantes mais próximos de O onde a solução explode, conclui-se que o 
problema de valor inical 





í d l = i-t + y 2 -ty 2 




[2/(0) = 


teria por solução 






y{t) = tan(í - -í 2 ) 


definida para 





í G]l - Vl + 7T,1 + VT 



<) 



(10) 



Suponha que temos uma equação diferencial da forma 

!-'(!)■ 

como, por exemplo, a equação Jr = sen(|). Estas equações dizem-se homogéneas. 
Como o segundo membro da equação depende apenas do quociente j, é natural 
fazer a substituição v = j, ou seja, y = wí. 
Mostre que esta substituição transforma ^ = /(|) na equação equivalente 






que é separável. 



Resolução: Supõe-se que t ^ 0. Se y(í) = íu(t), eníáo 



dy di> 

— = *~T + w 
aí aí 



dy _ f (y^ 



pelo que a equação -^ = / (| ) /Vca 






que é separável: 



1 



/(«)-« í 



para /(«) ^ t> 



D 



(11) 



Determine se cada uma das seguintes funções de t e y pode ser expressa como uma 
função de uma só variável |: 

(a) lny- lní + f_ 

(b) ^72+p- • 



t+y 
-y 



AMIV - FICHA SUPLEMENTAR 3 11 

Resolução: 

(a) Sim: fundindo os logaritmos e dividindo numerador e denominador da fracção por t, 
fica 

ln 2/ -lnt + ^ = lnf + i±|. 
t-y t 1 - f 

(b) Sim para t ^ 0: dividindo numerador e denominador da fracção por t s , fica 

y 3 + t 3 (|) 3 + i 

yt 2 + y 3 | + (|)3- 

Para t = 0, a função original vale 1 (ou tem limite 1 quando y — > 0), mas a substi- 
tuição para a variável | não é legal. 

D 

Comentário: Na alínea (b), a resposta estrita é "não". No entanto, pode ser útil saber 
exprimir a função em termos de uma só variável no domínio t ^ 0. <0 



Determine a solução geral da seguinte equação diferencial 



(12) 



,2 d v , . 2 
t M =2ty + y 



Sugestão: Exercício 10. 



Resolução: Para t ^ 0, divida-se a equação por t 2 e depois aplique-se o exercício 10: 

f = 2f + (f) 2 <=> tf t +v = 2v + v 2 



V _ 1 



í+í* - t parav + v 2 ^0 

J-^dv = fldt + c. 



Cálculo de uma primitiva de — p 

* 1)4-' 



v-\-v 2 



Decompõe-se , 2 em fracções simples: 



A i B 



1 = A + Av + Bv 



v(l-\-v) v l-\-v 

^^ A=\ and B = -1 



Logo, uma primitiva de } 2 é 



1 /" 1 

— dv — / dv = In \v\ — In |1 + v\ 

V J 1 + V 
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Continuando a resolução da equação: 

í —r^ dv = [ j dt + c <í=> In \v\ — In li + v\ = In líl 

J v\v z J t II I I II 



In 



= lnlíl 



l+v 

^^ tj- = kt onde k^O 

<í=> v — kt — ktv = O 
^> v(l - kt) = kt 

onde na última linha se substituiu v de volta por |. 

Apesar da técnica de resolução só se aplicar quando t ^ O, a extensão da função y acima 
para t = O é válida e permanece solução da equação. 

Falta estudar os casos em que v + v 2 = v(v + 1) = O (onde v = | ) excluídos no início 
da resolução: 

- quando v = O, obtém-se que y(t) = O, Vi, é solução, que pode ser descrita na forma 
acima para k = 0; 

- quando v = — 1, verifica-se, por substituição na equação diferencial, que y{t) = —t, 
Vi, também é solução: 

í 2 |(-í) = -í 2 = 2t(-t) + (-tf . 
Solução geral: 

kt 2 

y(t) = com k e M. ou y(t) = —t . 

1 — kt 

Intervalo de definição: Para k ^ O, o intervalo de definição de uma solução da primeira 
forma é 

] -oc,-[ ou ] 7 ,oo[ . 

k k 

Para k = O, o intervalo de definição de y(t) = O e R. O intervalo de definição da solução 

y(t) = —t também é R. 



Verificação: 



dy _ 2kt(l-kt)+k 2 t 2 

dt ~ (1-fct) 2 

,2dy_ _ 2kt 3 (l-kt)+k 2 t 4 

1 dt — (1-fcí) 2 



2ty + y 2 = 2kt3 ' kH4 



1-fcí -r (i-kt) 2 

_ 2kt 3 (l-kt) + k 2 t A _ , 2 dy _., 

~ (1-kt) 2 ~ l dt VK - 



D 



Comentário: A equação dada tem duas famílias de soluções definidas em subintervalos 
de R, cada uma parametrizada por k ^ 0, e tem ainda mais duas soluções definidas em 
todo o R. O 
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(13) 



Um objecto de massa m é lançado verticalmente a partir da superfície da Terra com 
uma velocidade inicial Vn. Considera-se um referencial em que o sentido positivo 
do eixo dos y's coincida com a direcção vertical apontando para cima, estando a 
origem sobre a superfície da Terra. Assumindo que não há resistência do ar, mas 
considerando a variação do campo gravitacional terrestre com a altitude, obtém-se 
a seguinte lei para a velocidade V(t) do objecto: 

dV mgR 2 

dt (y + R) 2 

onde R é o raio da Terra. 

(a) Seja V(t) = v(y(t)), onde v = v(y) é a velocidade como função da altitude y. 
Determine a equação diferencial satisfeita por v(y). 

(b) Calcule a chamada velocidade de escape, ou seja, calcule a menor velocidade 
inicial Vo para a qual o objecto não regressa à Terra. 

Sugestão: A velocidade de escape é determinada impondo que v(y) 
permaneça positivo. 



(14) 



Considere uma espécie com reprodução sexuada: cada membro da população ne- 
cessita de encontrar um parceiro para se reproduzir. Se N(t) for a população total 
no instante t, o número de encontros entre machos e fêmeas deve ser proporcional 
ao produto do número de machos pelo número de fêmeas. Como cada um destes 
números é proporcional a N(t), o número de nascimentos é proporcional a N 2 (t). 
Por outro lado, a taxa de mortalidade, é proporcional a N(t) pois não depende de 
encontros entre indivíduos. 

Conclui-se que a população N(t) satisfaz a equação diferencial 

dN , 

= bN 2 - aN , com a, b > 

dt 

Mostre que, se N(0) < f então N(t) — > quando t — > oo. Conclua que quando a 

população é inferior ao nível crítico | a população está em vias de extinção. 



Equações Exactas 



(15) 



Determine a solução geral da seguinte equação diferencial 



,'>, t , ir> ,. . , , o„.2A\dy 



2tsmy + y à e t + (r cosy + 3yV)-^ = 



Resolução: Esta equação não é linear nem separável. 
Teste de equação exacta: 
A equação 

2í sin y + y 3 e f + (t 2 cos y + 3y V) — = 
v v ' » „ ' dt 

M(t,y) N(t,y) 

é exacta se e só se existe uma função "potencial", F = F(t, y), tal que 

f = M(t,y) 



14 
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(por definição de equação exacta,). 

Uma vez que as expressões de M(t,y) e N(t,y) são continuamente diferenciáveis em 
todo o IR 2 , uma condição necessária e suficiente para a existência de uma tal função 
potencial F e 1 

dM dN 



Ora 



dy dt 

9M n „ 2 t dN 

— — = 2í cos y + 3y e = — - , 
oy ot 



pelo que se conclui que a equação é exacta. 

Resolução da EDO : 

Determina-se uma função "potencial" , F : 

^ = 2tsmy + y 3 e t 

ff = t 2 cos y + 3j/V 

F(t,y) = j(2tsmy + y 3 e t )dt + f(y) 

F(t,y) = J(t 2 cos y + 3y 2 e l ) dy + g(t) 

F(t,y) = t 2 smy + y 3 e t + f{y) 

F(t,y) = t 2 siny + y 3 e i + g(t) 

Escolhendo 

F(t,y)=t 2 smy + y 3 e t , 
a equação diferencial fica equivalente a 



±F(t,y(t)) = <=► F(t,y) 



Solução: 



í 2 sinj/ + j/V = c (*) 



t sin y + y ê = c . 



com c e 



Verificação: 



t 2 siny + y 3 ê 



e constante 



d (.2 



4 (í 2 sin í/ + 2/V) =0 



f (í 2 sinj/ + j/ 3 e < ) 



d_ (+2 

dy 



t siny + y e 



3„t\ dy _ 



dt 



2tsiny + y 3 e t + (t 2 cosy + 3y 2 e t )ifr = Q 



ok! 



D 



Comentário: Pelo teorema da função implícita, se for dada uma condição inicial, j/(ío) = 
yo. 3 expressão (*) determina y = y{t) em torno de t = to, desde que seja satisfeita a 
condição 

dF 

^-Oo,J/o) 7^ , 
dj/ 

ou seja, 

t 2 cosy + 3yy° ^0 . 



Cf. Análise Matemática III: condição para um campo vectorial (M(t, y), N(t, y)) ser o gradiante, \7F = 

F dF - 
dt ' dy ■ 



(llt ' §£)' de a| g uma função escalar F(í,y) 
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(16) 



Determine a 


solução 


do 


segu 


nte problema de valor inicial 

Í3í 2 + 4íy+(2y + 2í 2 )f = 

\y(o) = i • 






Resolução: A equação diferencial não é linear nem separável. 
Teste de equação exacta: 
A equação diferencial 

M 2 + Aty + (2y + 2t 2 ) ^ = 

M(t,y) N(t,y) 

é exacta se e só se existe uma função "potencial", F = F(t, y), tal que 

f = M(t,y) 

f = W,v) 

(por definição de equação exacta,). 

Uma vez que as expressões de M(t,y) e N(t,y) são continuamente diferenciáveis em 
todo o IR 2 , uma condição necessária e suficiente para a existência de uma tal função 
potencial F é 2 



Ora 



pelo que se conclui que a equação é exacta. 

Resolução da EDO : 

Para encontrar uma função "potencial" , F, resolve-se o sistema de equações 



dM 




dN 


dy 




~dt ' 


dM 
dy 


At 


dN 
~ ~dt 



dF 

dt 



dF 



= M 2 + Aty 



% = 2y + 2f 

F(t,y) = J(3t 2 +Aty)dt + f(y) 

F(t,y) = J(2y + 2t 2 )dy + g(t) 

F(t,y) = t 3 + 2t 2 y + f(y) 

F(t,y) = y 2 + 2t 2 y + g(t) 



Escolhendo 



F{t,y)=t 3 + y 2 + 2t 2 y 
a equação diferencial fica equivalente a 

d 



dt 



F(t,y(t)) = 



F(t,y) = c 



t 3 + y 2 + 2t 2 y = c . 



Cf. Análise Matemática III: condição para um campo vectorial (M(t, y), N(t, y)) ser o gradiante, VF = 

3F' 

dt ' dy- 



(llt ' §£)' de a| g uma função escalar F(t,y) 
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Condição inicial: 

A condição y(0) = 1 impõe que a constante real c satisfaça 

O 3 + l 2 + 2 • O 2 • 1 = c ^> c = 1 . 

A solução do problema de valor inicial é dada por 

t 3 + y 2 + 2t 2 y = 1 ^^ y 2 + 2t 2 y + í 3 - 1 = O 



«=► y(t) = -t 2 ± Vi 4 - í 3 + 1 . 

Para que y(0) = 1, escolhe-se o sinal + para a raiz quadrada. 
Intervalo de definição: 

- Para que y = y(t) esteja definida, o argumento da raiz quadrada tem que ser não- 
negativo. 

- Para que y = y{t) seja diferenciável, o argumento da raiz quadrada tem que ser positivo, 
o que acontece sempre: 

• quando t > 1, tem-se t 4 > t 3 , logo í 4 - í 3 + 1 > 1; 

• quando t < O, tem-se t 4 > O e -t 3 > O, logo í 4 - t 3 + 1 > 1; 

• quando O < t < 1, tem-se í 4 > O e i 3 < 1, logo í 4 - t 3 + 1 > 0. 
Conclui-se que o intervalo máximo de definição é IR. 

Solução do problema de valor inicial: 







y(t) = -t 2 

Verificação: 


+ Vi 4 - í 3 + 1 , para todo oíel. 


d u 
~dt 


o f | 4t 3 -3í 2 




- r ' 2Ví 4 -í a + l 


2y + 2t 2 


= 2Ví 4 -í 3 + l 


(2y + 2t 2 ) d J 


= -4íVí 4 - í 3 + 1 + 4í a - 3t 2 (•) 


3í 2 + Aty 


= 3í 2 - 4í a + 4íVí 4 - í 3 + 1 (**) 



(**) + (*) = O - ok! 



y(0) = -O + VO - O + 1 = 1 - ok! 



D 





Determine a constante real a para a qual a equação diferencial seguinte é exacta e 


(17) 


resolva-a: 

e at+y + 3t 2 y 2 + (2yt 3 + e at+y )^ = . 



Resolução: 

Determinação de a : 
A equação 



i at+y + 3t 2 y 2 + (2yt 3 + e at+y ) ^ = 
v ' » v ' dt 

M(t,y) N(t,y) 
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é exacta se e só se existe uma função "potencial", F = F(t, y), tal que 

f = M(t,y) 
f = N(t,y) 

(por definição de equação exacta/ 

Uma vez que as expressões de M(t,y) e N(t,y) são continuamente diferenciáveis em 
todo o M 2 , uma condição necessária e suficiente para a existência de uma tal função 
potencial F é 3 



dM 
dy 



dN 
~~dt 



Há então que escolher a de modo a satisfazer a equação 

e at+y + Qt 2 y = 6yt 2 + ^at+y ^^ Q = j _ 

Para que a equação seja exacta, a constante a deve ser 1. 

Resolução da equação para a = 1 : 

Quando a = 1, encontra-se uma função "potencial", F, resolvendo o sistema de equações 

% = é+y + zt 2 v 2 

f = 2yt 3 + e^ 

F(t,y) = f(é+y + 3t 2 y 2 )dt + f(y) 

F(t,y) = í{2yt z + e t +y)dy + g{t) 

F{t,y) = e t+ y+t 3 y 2 + f(y) 

F(t,y) = y 2 t 3 + e t+ y + g(t) 

Uma solução possível é 

F(t,y)=t 3 y 2 + e t+y . 
A equação dada é equivalente a 

d 



(lf 



F(t,y(t)) = 



*3„,2 



Solução: 



Verificação: 



tV 



t 3 y 2 



e t+y = c 



,t+y 



F(t,y) = c * 

(*) com c G 



í V + e^y = c 



t+y 



e constante 

f t {t 3 y 2 + é+y) = o 

í (*V + e t+ v) + £ (*V + e^y) § = 
3t 2 y 2 + e t+y + (2yt 3 + e t+ y)^- = -ok! 



D 



Cf. Análise Matemática III: condição para um campo vectorial (M(í, y), N(t, y)) ser o gradiante, VF = 
(^t"> lf )■ de a| g uma função escalar F(t, y). 
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Comentário: Pelo teorema da função implícita, se for dada uma condição inicial, y(to) = 
yo, a expressão (*) determina y = y(t) em torno de t = to, desde que seja satisfeita a 

condição 

dF 

-j-(to,yo)^0 , 
dy 

ou seja, 

2y íB + e to+yo + , 
É possível ver que para esta equação diferencial, as condições iniciais só são admissíveis se 
satisfizerem a condição do teorema da função implícita: 

Suponha-se (por redução ao absurdo) que existia uma solução com condição inicial 
(to, 2/o ) tal que 

2í/oío + e ío+yo = . 
A própria EDO avaliada em (ío,2/o) imporia também que 



e to+yo + 3t 2 y 2 = q _ 






A/esse caso, teria que serto^O e yo ^ 0. Então 






2y t 3 = 3t 2 y 2 <=^> 2t = 3y <= 


> yo = 


2 
= -t • 
3 


Voltando a substituir, ficaria 







0=2yotl + e t ^\ yo=¥o = ^4 + e^ , 

o que é impossível, pois a exponencial é sempre positiva e o primeiro termo nunca é 
negativo. <0> 



(18) 



Mostre que qualquer equação separável da forma 

dy_ = g{t)_ 
dt ~ f(y) 
é exacta. 



Resolução: Para f(y) ^ 0, 

dy g(t) . , , .dy 

-77 = JJ-T ^=> 9{t) -f{y) -77 = ■ 

dt f {y) ^-v^ > — „ — - dt 

M(t,y) N(t,y) 

A equação é exacta se e só se existe uma função "potencial", F = F(t,y), tal que 

f = M(t,y) 

f = W,V) 
ou seja, 

f = 9(t) 

% = -f(y) 

<=> F{t,y) = fg{t)dt-ff{y)dy + c 
Portanto, qualquer equação separável da forma 

dy_ = g(t) 
dt ~ f(y) 
é exacta; além disso, é equivalente, para f(y) ^ 0, a 

/ g(t) dt — f(y) dy = constante . 

D 
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Equações Redutíveis a Exactas 

(19) 



Resolva o seguinte problema de valor inicial 

Sty + y 2 + (t 2 + ty)y = , y(2) = 1 



Resolução: A equação não é linear nem separável. Também não é exacta porque: 

— {3ty + y 2 )=3t + 2y^2t + y = -{t 2 + ty) 

Assim, resta tentar encontrar um factor de integração p. A equação que p tem de verificar 
para ser um factor de integração é 

^-((3ty + y 2 )p) = -^((t 2 + ty)p) 

<=^ (3í + 2y)p + (3*1/ + y 2 )^- = (2í + y)p + (í 2 + ty)-£ 

Como não se aprende a resolver estas equações em AMIV, tenta-se achar um factor de 
integração que seja só função de y ou só função de t. Nesse caso a equação anterior 
simplifica-se porque uma das derivadas parciais de p é nula. 

A equação tem um factor de integração p = p{y) ? Em caso afirmativo, o factor de inte- 
gração tem de satisfazer a equação: 

SL(p(Zty + y 2 )) = ^{p{t 2 + ty)) 

^ ^(3*1/ + y 2 ) + p(y)(3t + 2y) = p{y){2t + y) 
dy 

éã t + y 



p 3ty + y 2 

Esta equação não tem soluções! Isto porque no termo esquerdo tem-se uma função de y 
e no termo direito uma função que depende de t e de y. Conclui-se que não existe um 
factor integrante que seja função só de y. 

A equação tem um factor de integração p = p(t) ? Em caso afirmativo, o factor de inte- 
gração tem de satisfazer a equação: 

— (p(3ty + y 2 )) = -{p(t 2 +ty)) (*) 

^ p(t)(3t + 2y) = p(t)(2t + y) + ^(t 2 + ty) 

í_ t + y 



p ty + t 2 

< > dt _ L 
p t 

Esta equação tem solução. Pode-se tomar por exemplo p{t) = t para t ^ 0. 

Resolução da equação usando o factor de integração: Multiplicando a equação pelo factor 

de integração obtem-se 

3t 2 y + y 2 t + (í 3 + t 2 y)y = (**) 
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que é equivalente ao problema que queremos resolver para t ^ 0. Como IR 2 é simplesmente 
conexo, a condição (*) é suficiente para a existência de uma função potencial para o campo 
vectorial (M, N) 

M(t,y) = 3t 2 y + y 2 t , N(t,y) = t 3 + t 2 y 

Sendo <p(t, y) a função potencial tem-se 

<P(t,y) = t 3 y + ±yH 2 + A(y) 

a ^ = t À + t 2 y { <P(t,y) = t 3 y + ±t 2 y 2 + B(t) 

Portanto, um potencial para o campo (M, N) é 




A equação (**) escreve-se 



(t>{t,y) = t 3 y+-t 2 y 2 



j t (Ht,y(t)) = o 



Portanto as soluções da equação verificam 

t 3 y + \t 2 y 2 = C 
e a constante pode ser determinada pelo valor inicial 

y {2) = 1 => 2 3 1 + -2 2 1 2 = C -^=> C = 10 
Ou seja, utilizando a fórmula resolvente, 

t 3 y + \t 2 y 2 -10 = 



-í 3 ± Vi 6 + 20í 2 

v = ê 



. , 20 

y=-t± ] /t2 + ¥ 

A condição inicial força o sinal "+" na equação anterior. 
Solução do problema de valor inicial: 



y{t) = -t+\Jt 2 + 2 j[ 

Esta é a solução do problema de valor inicial pretendido para í ^ 0. Como o limite de 
y{t) é infinito quando t — > 0, conclui-se que o intervalo máximo de definição da solução é 
]0, +oo[ (e portanto a condição que impusemos ( t ^ 0) é irrelevante para a solução deste 
problema de valor inicial). 
Verificação: 



20 
2/ (2) = -2+ A /4+ T = l 



Quanto à equação, fazendo u = Jt 2 + le notando que 



du_ 2t-f _ t 4 -20 
dt 2u t 3 u 
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tem-se 

.2 , 1.1 



3*1/ + y 2 + (t z + ty)y 

= 3t(-t + u) + (-t + u) 2 + (t 2 + t(-t + u)) ( -1 + 

= -3í 2 + 3íu + t 2 - 2tu + u 2 - tu + í 2 - -rr 

t z 

= - ofc! 



t 4 - 20 
t 3 u 



D 



(20) 



Resolva o seguinte problema de valor inicial 

y 3 + 2yt + (Ay 2 t + 2t 2 )y = , 2/(1) = 1 



A equação diferencial 



6y(í + y) + í(4í + 9y)^ = 



(21) 



admite um factor de integração da forma nity), ou seja, um factor fi que só depende 
do produto das variáveis ty. Determine-o e dê a solução da equação com y(l) = 1. 

Sugestão: A equação diferencial que dá fj, = fJ,(ty) pode ser escrita 
em termos de uma só variável v = ty. 



Resolução: 

Cálculo de um factor de integração: Uma função, u = n(ty), que nunca se anula é factor 
de integração da equação dada se e só se 



H(ty) 



6y(t + y) + *(4* + 9y) ^ 



= é equação exacta. 



Para tal, uma condição necessária é 
a 



% Wy)Qy{t + y)\ = & [u(ty)t(At + 9y)] 

tfi'(ty)Gy(t + y) + 6*/i(*</) + I2yfi{ty) = 

= yu'{ty)t{U + 9y) + /i{ty)8t + v{ty)9y 

(2t 2 y - 3ty 2 ) fi'(ty) + (3y - 2í) p(ty) = 



ty(2t-3y) -(2í-3y) 

Pondo v = ty, é suficiente resolver a seguinte equação separável com /i ^ e w ^ 0: 
vu>(v)-u(v) = ^ ^ = 1 

<í=> In |/u(w)| = In \v\ + c 

<^=> \lí{v)\ = k\v\ com k > 

<í=> fi(v) = kv com k ^ . 

Como basta um factor de integração, escolhemos k = 1, ou seja, 

fi(ty) = ty . 
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Resolução da EDO usando o factor de integração encontrado: Para tj^Oey^O, a 
equação dada é equivalente à seguinte equação exacta: 



ty 



6y(t + y) + t(U + 9?/)J 



= . 



Uma função "potencial", F(t,y), determina-se através de 

f = 6ty 2 (t + y) 

§ = t 2 y(4t + 9y) 

F(t,y) = J6ty 2 (t + y)dt + f(y) 

F{t,y) = ft 2 y{4t + 9y)dy + g{t) 

F(t,y) = 2íV + 3íV + f(y) 

F(t,y) = 2t 3 y 2 + 3t 2 y 3 +g(t) 

Uma solução possível é 

F(t,y) = 2t 3 y 2 + 3t 2 y 3 . 
A solução da equação diferencial é dada implicitamente por 

-F(t, y) = <^> F(t, y) = c ^ 2t 3 y 2 + 3t 2 y 3 = c . 

Verificação : Para ty ^ 0, 

2t 3 y 2 + 3t 2 y 3 é constante 

<=> é(2*V + 3tV)=0 

^ | (2t 3 y 2 + 3*V) + £ (2*V + 3*V) § = 



««=> íl/[6i/(í + »)]+ti/[í(4í + 9y)]^ = 

<«=► 6y(t + y)+ í(4í + 9y) f = - o/c/ 

Solução do problema de valor inicial: A condição y{\) = 1 impõe que a constante c satis- 
faça 

2-l-l + 3-l-l = c ^^ c = 5 . 
Pelo teorema da função implícita, como 

dy 
conclui-se que a solução do problema de valor inicial é dada por 

2t 3 y 2 + 3t 2 y 3 = 5 
para t numa vizinhança de to = 1 • D 



-(1,1) = t 2 y(4t + 9y)\ t=ly=l = 13^0 



Comentário: Esta solução poderia ser escrita explicitamente usando a fórmula resolvente 
para a equação do 3^- grau. (} 
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Existência, Unicidade e Extensão de Soluções 



(22) 



Determine uma solução contínua do problema de valor inicial 

í + y = 9(t) 

j/(0) = 

onde 

Í2 se < t < 1 

» í= in + Ci~ 

se t > 1 . 



Resolução: Para < í < 1, /x(í) = e* é um /actor de integração: 

^ K^y) = 2 e * 
^^ e*y = / 2e* dt + c 

^^ j/(í) = e" í (2e í + c) = 2 + ce _í 

Condição inicial : = y(0) = 2 + c => c = —2. 
Solução do problema de valor inicial para < t < 1: 

y(t) = 2 - 2e-* . 

Continuação da solução para í > 1: Quando t > 1, a equação é homogénea; para j/^0 
resolve-se pelo método de separação de variáveis: 

^^ f^dt = -fdt + c 

«=> !\dy = -t + c 

<í=> In |y| = — í + c 

<S=> |y(í)| = A;e _í onde fc > 

<S=> í/(í) = fce _í onde fc 7^ . 

/\ função y{t) = 0, Vi, também é solução. Logo, a solução geral da equação diferencial 
para t > 1 é 

y(í) = fce _í onde fcel. 

Para obter uma solução contínua do problema posto, a solução para t > 1 deve satisfazer 
uma condição inicial em t = 1 que a faça coincidir com o valor da solução no intervalo 
< í < 1: 

y(l) = 2-2e- 1 . 

Assim impõe-se ke~ x = 2 — 2e _1 , ou seja, k = 2e — 2. 
So/tvção : 

_ 12-26"* , < í < 1 

V{t) ~ \ (2e - 2)e-* , í > 1 . 



24 
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Verificação: 



dy_ 
dt 



2j[2-2e - *] = 2e-* , O < í < 1 

^[(2e-2)e-*] = (2-2e)e-* , í>l. 



dt ^ y 



he"* + 2- 2e _i = 2 , O < í < 1 
\(2-2e)e _f + (2e-2)e _f = , t > 1 
= g(í) para í ^ 1 - o/c/ 



D 



Comentário: Esta "solução" contínua não é de classe C 1 , porque os limites laterais de 



J em í = 1 são diferentes: 



lim § = lim 2e _í = 2e _1 
t->i- t-»i- 



dy _ 



lim ^ = lim (2 - 2e)e _ * = (2 - 2e)e" 1 . 
í—i+ dt í—i+ 



Portanto, esta "solução" não é solução do problema de valor inicial no sentido estrito da 
definição de "solução de equação diferencial" onde se impõe ter primeira derivada contínua 
(i.e., ser de classe C 1 ). (} 



Prove que y(t) = — 1, Vi G 



(23) 



é a única solução do problema de valor inicial 

f dy _ 

) dt 

1 2/(0) = 



dt *(i + y) 



■i . 



Resolução: /A função y(t) = — 1, Vt G R, é claramente solução do problema dado porque 

d 



dt 



(-1) = e í(l + (-l)) = 0, VígR. 



Como f(t,y) = í(l + j/) é uma função de classe C 1 definida em R 2 , ela é localmente 
lipschitziana relativamente a y em M 2 . Pelo teorema de Picard, qualquer problema de 
valor inicial 

ÍJ = f(t,y) = t(l + y) 

\y(to) = yo 

tem uma única solução. Em particular, y(t) = — 1, Vi G R, é a única solução de 

/* = *(! + ») 
[2/(0) = -1 • 

D 



(24) 



Determine 


todas 


as 


soluções 


de 


classe 
dy _ 
dt ~ 


C L do 

y/y > 


problema 
1/(0) = 


de 

= 


valor 


inicial 
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Mostre que existe uma solução de classe C 1 para o problema de valor inicial 



(25) 



2/(0) = -1 
diferente da solução y(í) = — 1, Vi 6 M. 

Sugestão: Para |y| < 1 a equação pode ser resolvida como equação 
separável. Obtenha dessa maneira um par de soluções com a propri- 
edade y(t) — > —1 quando í — > 0+ ou í — > _ , cole-as e estenda ao 
resto de R como constante. 

Explique porque é que isto não contradiz o teorema de Picard. 



Resolução: De acordo com a sugestão, resolve-se a equação diferencial para \y\ < 1 
como equação separável: 



d J = tVT^ <=> -£_ = * 






arcsin y = -^ + c 



y(t) = sin(Ç + c) . 



Para que y(0) = — 1, escolhemos c = — |. 

Para que |y| < 1, o argumento do seno tem que ser 



t' 2 



-f < y- f < f ^^ 0<í"<2^ 



te] - V27T,0[ OU* e]0,V27r 



Sejam 



y-(í)=sin(Ç-f), íe]-v^F,0[, 

y+(í) = sin(Ç - f ) , ÍG]0,v^F[. 



Estas soluções podem-se colar em t = e podem-se estender para t > V27r ou t < 
\Fhi como sendo a constante 1. 



Solução: 



() <M fsin(V-f) , |í|<V2^ 

|1 , 1*1 > v^ 
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Note-se que esta função é continuamente diferenciável em todo o R (em particular, a 
derivada em t = — v / 2~7r ou em t = O ou em t = VÕjt é 0). 



Verificação: 



ícos(Ç- f) , \t\ < V^TT 



dy \ í-CUí^y ~~ 2> > 

dt [O, \t\ > y/2ir 



ty/l - sin 2 (f - f) , |í| < v^2tt 

O , \t\> \Ph\ 



t^T^p = rv— "---' 



-ok! 



ícos(Ç- |) , \t\ < y/^r 



[O , \t\ > V2tt 

Relação com o teorema de Picard: A equação diferencial é da forma 

djL = f{t,y) 
dt Jy,y ' 



com f(t,y) = t\/l — y 2 definida para t el e y e [—1, 1]. 

Pelo teorema de Picard, a solução do problema de valor inicial 

\y{to) = yo 

existe e é única se f for localmente lipschitziana em relação a y num domínio contendo 
(*o,2/o)- 

A/este caso, a função f(t,y) = t\/l — y 2 não é localmente lipschitziana em relação a 
y para qualquer domínio contendo (ío,2/o) = (0,-1)- A condição falha exactamente em 
y = — 1 porque 

lim — — = lim — - = oo para t ph , t ^ , 

y— i dy y— i yi^2 

o que implica que 

lim t = oo para í»0 ,íf O . 

«— 1 y-(-l) 

/.ogo, qualquer que seja o intervalo [—1,-1 + e], não pode existir uma constante L £ 
satisfazendo 4 

\f(t, y) - f(t, -1)| < L £ \y + 1| , Vy e [-1, -1 + e] , 

para t « 0, í ^ 0. D 



Comentário: 

• /4 função f(t,y) é localmente lipschitziana em relação a y para y e] — 1, 1[, porque 
a derivada parcial -^- é contínua para y e] — 1, 1[. 

• /4 função y{t) = sin(y — |), Vi e IR, não é solução da equação diferencial. De 

facto, a igualdade vi — sin 2 x = cos a; usada na verificação só é verdadeira quando 
— § + 2fc7r < x < § + 2fc7r para algum fceZ. Para outros valores de x, o coseno é 



negativo e tem-se yl — sin 2 a; = — cos x. 



Cf. definição de função localmente lipschitziana. 



(26) 
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Mostre 


que a 


solução do 


segu 


inte 


pro 


blema 


de 


valor inicial 


existe, 


é 


única 


e 


está 


definida 


para 


< t < 1: 




y = 


-y 2 


+ cos 


t\ 


y(o) = 


= 













Resolução: Esta equação não pode ser resolvida explicitamente pelos métodos estudados. 
No entanto, a resolução explícita não é necessária para responder à questão. 

A função f(t,y) = y 2 + cosi 2 é de classe C 1 em M 2 e portanto, o teorema de Picard 
garante a existência e unicidade da solução do problema de valor inicial. O teorema garante 
também que a solução do problema de valor inicial pode ser prolongada a um intervalo 
máximo de definição ]a,b[ (contendo 0) tal que quando t —> a + ou t —> b~ , (t,y(t)) 
tende para a fronteira do domínio de f. Uma vez que f está definida em todo o R 2 , 
isto significa que (t,y(t)) —> oo nos extremos do intervalo de definição. Em particular, se 
a (respectivamente b) for finito então a solução explode para t = a (respectivamente 
t = b), isto é, y(t) —> oo quando t —> a (respectivamente t —> b). 

Assim, para mostrar que o intervalo de definição da solução contém [0, 1] é suficiente 
mostrar que a solução y{t) não explode para t < 1. Uma maneira de fazer isto é arranjar 
dois problemas de valor inicial 

^ = /(*»«) - «(0) = 

dv , 

— = g(t,v) , u(0) = 

cujos intervalos de definição contenham [0, 1] e que verifiquem 

u{t) < y{t) < v{t) 

Para que esta última condição seja verificada basta 5 que 

f(t,y) <y 2 + cos t 2 <g(t,y) 

Uma vez que, para < t < 1 se tem 

< y 2 + cos t 2 < y 2 + 1 

pode-se considerar os problemas de valor inicial 

|.0, ,(0)=0 

f = „s + l,„«» = 

O primeiro tem solução constante u(t) = 0. Quanto ao segundo, trata-se de uma equação 
separável 

1 dv 



v 2 + l dt 
<í=> — (arctanu) = 1 
<í=> v(t) = tan(í + c) 
em que c é uma constante real. A condição inicial implica c = 0, portanto a solução é 

v(t) = taní 



5 



Ver Proposição 3.2.11, página 160, do livro "Equações Diferenciais Ordinárias" por Fernando 



Pestana da Costa. 



28 
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com intervalo máximo de definição ] — §,§[■ Conclui-se que 

O < y(t) < taní, para O < t < 1 
o que mostra que o intervalo de definição de y{t) contém [O, 1]. 



D 



(27) 



Mostre que as soluções dos seguintes problemas de valor inicial existem, são únicas 
e estão definidas nos intervalos indicados: 

(a) y = y + e~ y + e"* , y(0) = , < t < 1 

(b) y = e"* + ln(l + y 2 ) , y(0) = , < t < +cx) 



Campos de Direcções 



(28) 



Esboce o campo de direcções e trace os respectivos tipos de solução da equação 
diferencial 

-r- = y{y-2) . 

dt yvy J 



Resolução: Para cada c e IR, o conjunto dos pontos (t, y) e M? onde o gráfico da solução 
y{t) tem declive -^ = c, é determinado pela equação 

y(y - 2) = c . 



Casos especiais: 

c< -1 
c= -1 



c = - 

c = 

c= 1 

c = 3 

c = 8 



y 2 - 2y - c = 
y 2 - 2y + 1 = 

y 2 - 2í/ + \ = 
y = ou y = 2 



V 
V 2 



2y - 1 = 
2y - 3 = 



y z - 2y - 8 = 



y = 1 ± vT+c é impossível 

y = i 

2/ = l±V2 
y = l±2 
y = l±3 
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Esboço do campo de direcções: 
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Traçado dos tipos de solução: 




D 



(29) 



Esboce os campos de direcções e trace os vários tipos de solução para as seguintes 
equações diferenciais 

(a) y = y{y 2 - 1) 

(b) y = t 2 +y 2 



